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FIELD OF THE INVENTION 

The present invention relates to fluid process 
control systems. In particular, the present invention 
relates to diagnostic systems for fluid flow in 
process control systems. 

BACKGROUND OF THE INVENTION 

Fluid flow meters are used in industrial process 
control environments to measure fluid flow and provide 
flow signals for flow indicators and controllers. 
Inferential flow meters measure fluid flow in a pipe 
by measuring a pressure drop near a discontinuity 
within the pipe. The discontinuity (primary element) 
can be an orifice, a nozzle, a venturi, a pitot tube, 
a vortex shedding bar, a target or even a simple bend 
in the pipe. Flow around the discontinuity causes both 
a pressure drop and increased turbulence. The pressure 
drop is sensed by a pressure transmitter (secondary 
element) placed outside the pipe and connected by 
impulse lines or impulse passageways to the fluid in 
the pipe. Reliability depends on maintaining a correct 



calibration. Erosion or buildup of solids on the 
primary element can change the calibration. Impulse 
lines can become plugged over time, which also 
adversely affects calibration. 

Disassembly and inspection of the impulse lines 
is one method used to detect and correct plugging of 
lines. Another known method for detecting plugging is 
to periodically add a "check pulse" to the measurement 
signal from a pressure transmitter. This check pulse 
causes a control system connected to the transmitter 
to disturb the flow. If the pressure transmitter fails 
to accurately sense the flow disturbance, an alarm 
signal is generated indicating line plugging. Another 
known method for detecting plugging is sensing of both 
static and differential pressures. If there is 
inadequate correlation between oscillations in the 
static and differential pressures, then an alarm 
signal is generated indicating line plugging. Still 
another known method for detecting line plugging is to 
sense static pressures and pass them through high pass 
and low pass filters. Noise signals obtained from the 
filters are compared to a threshold, and if variance 
in the noise is less than the threshold, then an alarm 
signal indicates that the line is blocked. 

These known methods rely on providing static 
pressure sensors or disassembly of the flow meter, 
increasing complexity and reducing reliability. These 
known methods do not provide for diagnosing the 
condition of the primary element. There is thus a need 
for a better diagnostic technology providing more 
predictive, less reactive maintenance for reducing 
cost or improving reliability. 

SUMMARY OF THE INVENTION 



A flow diagnostic system is disclosed for 
coupling to a primary flow sensing element via impulse 
lines. The flow diagnostic system may include a 
pressure transmitter that couples to the impulse lines 
and generates digital pressure data representing 
pressure . 

A control system receives the pressure data and 
provides the pressure data and real time clock 
readings associated with pressure data to a diagnostic 
application stored in the flow diagnostic system. 

The diagnostic application includes a first 
algorithm that calculates a difference between the 
pressure data and a moving average of the pressure 
data. 

The diagnostic application also includes a second 
algorithm that receives the difference and calculates 
a trained data set of historical pressure data during 
a training mode and calculates a current pressure data 
set during a monitoring mode and then generates 
diagnostic data as a function of the current pressure 
data set relative to the historical pressure data, 
which indicates changes in the condition of flow 
sensing. 

The diagnostic application includes a third 
algorithm that generates a report indicating the 
diagnostic data. 

The diagnostic application can be provided over a 
network from an application service provider (ASP) or 
other source and can run on the control system, the 
ASP or another remote computer. 



BRIEF DESCRIPTION OF THE DRAWINGS 



FIG. 1 is a schematic illustration of a flow 
diagnostic system that diagnoses the condition of 
impulse lines and/or a primary flow element. 

FIG. 2 is a schematic illustration of a first 
embodiment of a control system that is part of a flow 
diagnostic system. 

FIG. 3 is a schematic illustration of a second 
embodiment of a control system that is part of a flow 
diagnostic system. 

FIG. 4 is a block diagram of a diagnostic 
application . 

FIG. 5 is a block diagram of a diagnostic 
application that diagnoses the condition of impulse 
lines . 

FIG. 6 is a block diagram of a diagnostic 
application that diagnoses the condition of a primary 
flow element . 

FIG. 7 is a flow chart of a process diagnosing 
the condition of impulse lines. 

FIG. 8 is a flow chart of a process diagnosing 
the condition of a primary element. 

FIG. 9 is a flow chart of a process of diagnosing 
the condition of both impulse lines and a primary 
element . 

FIG. 10 is a graph of amplitude versus frequency 
versus time of a process variable signal. 

FIG. 11 is a block diagram of a discrete wavelet 
transformation . 

FIG. 12 is a graph showing signals output from a 
discrete wavelet transformation. 

FIG. 13 is a diagram showing a simplified neural 
network . 



FIG. 14 is a diagram showing a neural network 
used to provide a residual lifetime estimate. 

FIG. 15 is a graph of residual life versus time. 

FIG. 16 illustrates a diagnostic fluid flow meter 
that has a pitot tube for a primary element. 

FIG. 17 illustrates a block diagram of a computer 
platform communicating with a transmitter via a field 
network . 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

In the present invention, a flow diagnostic 
system is provided . that can diagnose the condition of 
either the primary element or impulse lines connected 
to a pressure transmitter. A diagnostic application is 
downloadable over a network from an application 
service provider (ASP) , or can be obtained from a 
computer-readable medium such as a CD-ROM or removable 
disc. The diagnostic application can run on the 
control system, a remote computer or the ASP and 
provide a diagnostic report. The diagnostic 
application runs on a processor that is high powered 
relative to the low power embedded processor found in 
the pressure transmitter. With the use of the high 
powered processor, sophisticated diagnostics can be 
performed in real time and provide prompt reports to a 
plant operator about the condition of primary elements 
or impulse lines or both. There is no need to use 
processing time on the processor imbedded in the 
transmitter for diagnostics, and the transmitter can 
perform its measurement tasks rapidly. 

FIG. 1 is a schematic illustration of a 
generalized example of a flow diagnostic system 100 
that diagnoses the condition of impulse lines 104 
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and/or a primary flow element 106 placed in a fluid 

piping system 108. The impulse lines 104 and the 

primary element 106 are referred to collectively as a 

"pressure generator . " 

The term "pressure generator" as used in this 

application means a primary element (e.g., an/orifice 

plate, a pitot tube, a nozzle, a venturi , >a shedding 

bar, a bend in a pipe or other flow .discontinuity 

adapted to cause a pressure drop in/rlow) together 

with impulse pipes or impulse passageways that couple 

the pressure drop from locations'^ near the primary 

element to a location outside/ the flow pipe. The 

spectral and statistical characteristics of this 

/ 

pressure presented by ^rthis defined "pressure 

generator" at a locationyOutside the flow pipe to a 

connected pressure transmitter 82 can be affected by 

the condition of the primary element as well as on the 

condition of the impulse pipes. The connected pressure 

transmitter can beya self-contained unit, or it can be 

fitted with remote seals as needed to fit the 

application. A/flange on the pressure transmitter 102 

(or its remote seals) couples to a flange adapter on 

the impulse lines 104 to complete the pressure 

connections in a conventional manner. The pressure 
/ 

transmitter 102 couples to a primary flow element 106 

via impulse lines 104 to sense flow. Primary element 



106 
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as illustrated, is an orifice plate clamped 



between pi pe flanges 105. ^_ 

*7 Pressure transmitter 102 is coupled to the 
30 impulse lines 104 and provides pressure data along 
line 110 to a control system 112. Pressure transmitter 
102 is preferably a differential pressure transmitter. 
Line 110 is typically a conventional 4-20 mA loop that 
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provides all of the electrical energization for 
transmitter 102 and also carries digital pressure data 
representing the pressure sensed by pressure 
transmitter 102. Pressure transmitter 102 includes an 
5 embedded microprocessor that generates the digital 
pressure data, however, the computing power of the 
embedded microprocessor is limited by the power 
available from the 4-20 mA loop. The digital pressure 
data is typically transmitted in a standard industry 
10 format using a digital protocol such as HART, 
Foundation Fieldbus, Profibus, CAN or the like. 

Control system 112 receives the digital pressure 
data on line 110 and providing the pressure data and 

real time clock readings associated with pressure data 

Lip 

§k.s 1 15 to a diagnostic application stored in the flow 

diagnostic system 100. 

As explained in more detail below, the diagnostic 

ru 

«, application can be a application 134 that is stored in 

0 and running on a processor in control system 112, an 

In 

20 application 130 stored in and running on a processor 
at an application service provider 126, an application 
132 stored in and running on a processor in a remote 
computer 122 or any combination of these locations. In 
each of these locations, relatively high power 
25 computing is available in comparison to the limited 
computing power available from the embedded processor 
in transmitter 102. An arrangement where the 
diagnostic application is stored on one computer (112, 
122 or 126) and running on another one of the 
30 computers (112, 122, 126) is also contemplated. 

The diagnostic application (134, 130, 132 or a 
combination of 134, 130 or 132) includes a first 
algorithm calculating a difference between the 
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pressure data and a moving average of the pressure 
data. The diagnostic application further includes a 
second algorithm receiving the difference and 
calculating a trained data set of historical pressure 
data during a training mode and calculating a current 
pressure data set during a monitoring mode and 
generating diagnostic data as a function of the 
current pressure data set relative to the historical 
pressure data indicating changes in the condition of 
flow sensing. The diagnostic application also includes 
a third algorithm generating a report indicating the 
diagnostic data. The report can be presented to an 
user at control system monitor 114 , a printer 124 or 
monitor associated with remote computer 122 or a 
report 128 generated by the application service 
provider 126. It is appreciated that the report can 
take a variety of forms or displays that communicate 
the diagnostic data to the user. 

The control system 112 couples along line 118 to 
a network 116. Network 116 also couples along line 120 
to the remote computer 122 , and couples along line 127 
to the application service provider 126. Typically, 
the remote computer 122 is located in the offices of 
the plant operator, while the line 127 is a long 
distance connection, typically over the internet, to 
an application service provider 126 at any desired 
location. 

The network 116 can provide the diagnostic 
application to the control system 112 or the remote 
computer 122. The network can also provide the 
pressure data and associated real time clock readings 
to the diagnostic application running at the remote 
computer 122 or the application service provider 126. 



The illustration in FIG. 1 is a generalized 
illustration, and not all diagnostics systems include 
all of the features shown in FIG. 1. 

FIG. 2 is a schematic illustration of a first 
embodiment of a control system 2 06 that is part of a 
flow diagnostic system such as flow diagnostic system 
100 in FIG. 1. 

Pressure transmitter 200 is coupled via impulse 
lines 201 to an integral orifice plate 202 in a fluid 
flow system. Pressure transmitter 200 couples to 
control system 206 along line 204. The pressure 
transmitter generates digital pressure data (P) and 
associated digital real time clock readings (RTC) as 
illustrated at 205. The real time clock readings 
provided by the transmitter 200 with each pressure 
reading indicate the relative time that each pressure 
reading is taken by the transmitter 200. In an 
instance where the transmitter generates pressure 
readings at a generally constant rate, the real time 
clock reading can be as simple as a sequential number 
for each reading. In an instance where the pressure 
reading are somewhat more irregularly spaced in time, 
the real time clock reading can be an approximate time 
that each reading taken. 

The control system 206 receives the digital 
information 205 at a transmitter bus interface 208. 
Transmitter bus interface 208 formats the data 
properly for communication on the bus 210 that is 
internal to the control system 206. The information 
205 is used by process control software 214 to couple 
control signals along line 215 to a valve or other 
control device. Control system 206 includes a network 
interface 212 that can pass the digital information 
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205 along to a network on line 218 in a case where the 
diagnostic application is not running on the control 
system 206. In a case where the diagnostic application 
is running on the control system 206, the network 
5 interface 212 can receive the diagnostic application 
from line 218 and format it for transmission over the 
bus 210. The network interface 212 is bidirectional. 
The control system 206 includes a processor, random 
access memory, non-volatile memory and I/O devices. 
10 In FIG. 2, the transmitter 200 provides real time 

clock readings. As explained below in connection with 
FIG. 3, a flow diagnostic system can also operate with 
a transmitter that does not provide real time clock 
readings . 

15 FIG. 3 is a schematic illustration of a second 

embodiment of a control system 3 06 that is part of a 
flow diagnostic system such as flow diagnostic system 
100 in FIG. 1. 

Pressure transmitter 300 is coupled via impulse 

20 lines 3 01 to an integral orifice plate 3 02 in a fluid 
flow system. Pressure transmitter 300 couples to 
control system 306 along line 304. The pressure 
transmitter generates digital pressure data (P) as 
illustrated at 305. Transmitter 300 may or may not 

25 generate real time clock readings, depending on the 
design of the transmitter. 

In an instance where the transmitter 300 
generates real time clock readings, these reading are 
coupled through transmitter bus interface 3 08 to real 

30 time clock circuit 318 via a bus 310 that is internal 
to control system 306. Real time clock circuit 318 
then generates a corresponding synchronized real time 
clock reading that is synchronized with other real 



time clock readings in the control system 306. This 
synchronization allows for synchronized comparison of 
data from multiple transmitters, for example, in a 
case where differential pressure is measured using two 
absolute pressure transmitters. 

In an instance where transmitter 3 00 does not 
generate any usable real time clock readings, then 
real time clock circuit 318 generates synchronized 
real time clock readings based on the time that the 
each pressure data point 305 is received by real time 
clock circuit 318. 

In either instance , the pressure data and 
synchronized real time clock readings can be used by 
process control software 314 to couple control signals 
along line 315 to a valve or other control device. 
Control system 306 includes a network interface 312 
that can pass digital information (pressure data 305 
and synchronized real time clock readings) along to a 
network on line 318 in a case where the diagnostic 
application is not running on the control system 306. 
If the network connected to line 318 is busy or 
otherwise unavailable, then the digital information 
can be temporarily stored in storage device 320 and 
later transmitted via the network to a diagnostic 
application running on a remote computer. This storage 
320 can also be used in instances where it is desired 
to use batch processing for diagnostics. 

In a case where the diagnostic application is 
running on the control system 306, the network 
interface 312 can receive the diagnostic application 
from line 318 and format it for transmission over the 
bus 310 internal to control system 306. 



FIG. 4 is a block diagram of a diagnostic 
application 432 running on a remote computer 422. 
Pressure data and real time clock reading from a 
pressure transmitter are provided on line 320 from a 
network 316. 

A first algorithm 440 stored in computer 422 
calculates a difference between the series of digital 
pressure data points received from network 316 and a 
moving average of the series of digital pressure data 
points. Real time clock readings associated with each 
digital pressure data point are used to ensure that 
the digital pressure data points, which may be delayed 
or out of order after passing through network 316, are 
reassembled in a correct time sequence at the input to 
algorithm 440. The difference calculated in algorithm 
44 0 is coupled along line 44 2 to a second algorithm 
444 . 

The second algorithm 444 is also stored in the 
computer 422. The second algorithm 444 receives the 
difference calculated by algorithm 440 and calculates 
a trained data set of historical data during a 
training mode and calculates a current data set during 
a monitoring mode and generates diagnostic data 446 as 
a function of the current data set relative to the 
historical data indicating changes in the condition of 
the pressure generator (not illustrated in FIG. 4) . 

A third algorithm 448 receives the diagnostic 
data 446 and generates a report 450 indicating the 
diagnostic data to a user. The report 450, as 
explained above, can take forms such as a monitor 
display or a printed report . 
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FIG. 5 is a block diagram of a diagnostic 
application running on a remote computer 522 that 
diagnoses the condition of impulse lines. 

In FIG. 5, a difference (subtract moving average) 
algorithm 540 receives digital pressure data from 
network 516 via line 520 and generates difference data 
on line 542 representing a sensed pressure minus a 
moving average of the sensed pressure. A calculate 
algorithm 552 receives the difference data from line 
542 and calculates a trained output 553 of historical 
data obtained during a training mode or time interval. 
After training, calculate algorithm 554 calculates a 
monitor output 555 of current data obtained during a 
monitoring mode or normal operation time of the fluid 
flow meter. 

. In FIG. 5 , a diagnostic algorithm 556 receives 
the trained output 553 and the monitor output 555 and 
generates a diagnostic data output 558 indicating a 
current condition of the pressure generator relative 
to an historical condition. In FIG. 5, calculate 
algorithm 554 stores the historical data in memory for 
later use in the comparison. 

In difference algorithm 540, the moving average 
is calculated according to the series in Eq. 1: 

m 

Aj = S ( p j+k) (W k ) Eq. 1 

k=0 

where A is the moving average, P is a series of 
sequentially sensed pressure values, and W is a 
numerical weight for a sensed pressure value, m is a 
number of previous sensed pressure values in the 
series. Provision can also be made in difference 
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circuit 54 0 to filter out spikes and other anomalies 
present in the sensed pressure. In FIG. 5, the 
historical data comprises statistical data, for 

example, the mean (ji) and standard deviation (<y) of 
the difference output or other statistical 
measurements, and the diagnostic output 558 indicates 
impulse line plugging. The diagnostic application 
switches between a training mode when it is installed 
and a monitoring mode when it is in use measuring flow 
as illustrated by switch 550. The calculate algorithm 
554 stores historical data in the training mode. The 
diagnostic data output 558 indicates a real time 
condition of the pressure generator. 

In FIG. 5 , statistical data, such as the mean ^ 

and standard deviation a, are calculated based on a 
relatively large number of data points or flow 
measurements. The corresponding sample statistical 
data, such as sample mean X and sample standard 
deviation s, are calculated from a relatively smaller 
number of data points. Typically, hundreds of data 
points are used to calculate statistical data such as 

jx and a/ while only about 10 data points are used to 
calculate sample statistical data such as X and s. The 
number of data points during monitoring is kept 
smaller in order to provide diagnostics that is real 
time, or completed in about 1 second. Diagnostic 
algorithm 556 indicates line plugging if the sample 
standard deviation s deviates from the standard 

deviation g by a preset amount, for example 10%. 

FIG. 6 is a block diagram of a diagnostic 
application running on a remote computer 622 that 
diagnoses the condition of a primary flow element. The 
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diagnostic application of FIG. 6 is arranged generally 
along the lines of the diagnostic application of FIG. 
5. A difference algorithm 640 receives digital 
pressure data from network 616 via line 620. The 
diagnostic application switches between a training 
mode and a monitoring mode by switch 650. 

In FIG. 6, the diagnostic output 658 indicates a 
condition of the primary element, while in FIG. 5, the 
diagnostic output 558 indicates a condition of the 
impulse lines. In FIG. 6, calculate algorithms 652, 
654 calculate and store data on power spectral density 
(PSD) of the difference output on line 642 rather than 
statistical data which is used in FIG. 5. 

The power spectral density data is preferably in 
the range of 0 to 100 Hertz. The center frequency of a 
bandpass filter can be swept across a selected range 
of frequencies to generate a continuous or quasi- 
continuous power spectral density as a function of 
frequency in a manner that is well known. Various 
known Fourier transforms can be used. 

Power spectral density, Fi, can also be 
calculated using Welch's method of averaged 
periodograms for a given data set . The method uses a 
measurement sequence x(n) sampled at fs samples per 
second, where n = 1, 2, . . .N . A front end filter with 
a filter frequency less than fs/2 is used to reduce 
aliasing in the spectral calculations. The data set is 
divided into Fjc,i as shown in Eq. 2: 



M 



Fk.i = (i/M) IE x k (n)e" j2niAfn 



Eq.2 
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There are Fk,i overlapping data segments and for 
each segment, a periodogram is calculated where M is 
the number of points in the current segment . After 
all periodograms for all segments are evaluated, all 
of them are averaged to calculate the power spectrum: 



Once a power spectrum is obtained for a training 
mode, this sequence is stored in memory, preferably 
EEPROM, as the baseline power spectrum for comparison 
to real time power spectrums. Fi is thus the power 
spectrum sequence and i goes from 1 to N which is the 
total number of points in the original data sequence. 
N, usually a power of 2, also sets the frequency 
resolution of the spectrum estimation. Therefore, Fi 
is also known as the signal strength at the i th 
frequency. The power spectrum typically includes a 
large number points at predefined frequency intervals, 
defining a shape of the spectral power distribution as 
a function of frequency. 

In the detection of the primary element 
degradation, a relatively larger sample of the 
spectral density at baseline historical conditions and 
a relatively smaller sample of the spectral density at 
monitoring conditions are compared. The relatively 
smaller sample allows for a real time indication of 
problems in about 1 second. An increase in the related 
frequency components of the power spectrum can 
indicate the degradation of the primary element. Using 



Fi = (1/L) 




Eq. 3 



k=l 



orifice plates as primary elements, for example, 
changes as high as 10% are observed in spectral 
components when the orifice plate is degraded to a 
predetermined level. The amount of change can be 
adjusted as needed, depending on the tolerable amount 
of degradation and the type of primary element in use. 
The amount of change needed to indicate a problem is 
arrived at experimentally for each type of primary 
element arrangement. Fuzzy logic can also be used to 
compare the many points of the power spectrums . 

The pressure data is typically a calibrated 
output, and the diagnostic data indicates if the 
pressure generator is out of calibration. 

FIG. 7 is a detailed flow chart of a process 
diagnosing the condition of impulse lines. 

The algorithm starts at 702. A moving average is 
subtracted from differential pressure data as shown at 
704 to calculate a difference. During a train mode, 
historical data on the calculated difference is 
acquired and stored at 706 as statistical data ^ and 
or, for example. During an operational MONITOR mode, 
current data on the difference is acquired and stored 
at 708 as statistical data X and s. The smaller sample 
of current data is compared to the larger sample of 
the historical data to diagnose the condition of the 
impulse lines. Comparisons 710 of historical and 
current statistical data are made at 714, 716, 718 and 
a selected diagnostic output is generated at 730, 732, 
734 as a function of the comparisons made at 712, 714, 
716, 718 respectively. 

After completion of any diagnostic output, the 
process loops back at 720, 722, 724, 726 or 728 to 
repeat the monitor mode diagnostics, or the 
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transmitter can be shut down until maintenance is 
performed. If the diagnostic process itself fails, an 
error indication is provided on the diagnostic output 
at 736. In the method of diagnosis illustrated in FIG. 
7, the historical data set comprises statistical data 
such as data on the mean (ji) and standard deviation 
(a) of the calculated difference; the current data set 
comprises current sample statistical data, such as the 
sample average (X) and sample deviation (s) of the 
calculated difference. The sample deviation (s) is 
compared to the standard deviation (a) to diagnose 
impulse line plugging, for example. Other known 
statistical measures of uncertainty, or statistical 
measures developed experimentally to fit this 
application can also be used besides mean and standard 
deviation. When there is an unusual flow condition 
where X_is much different than jj,, the diagnostics can 
be temporarily suspended as shown at 712 until usual 
flow conditions are reestablished. This helps to 
prevent false alarm indications. 

In FIGS. 1-7, the transmitter generates a 
calibrated output and the flow diagnostic system 
generates a diagnostic output that indicates if the 
pressure generator is out of calibration. In FIGS. 1- 
7, the primary element can comprise, for example, a 
simple pitot tube or an averaging pitot tube as 
illustrated below in the example shown in FIG. 16. The 
primary element and impulse pipes can be combined in 
an integral orifice such as illustrated in FIGS 2-3. 
An orifice plate adapted for clamping between pipe 
flanges is illustrated in FIG. 1. The primary element 
can comprise a venturi, a nozzle or an orifice (not 
illustrated) . A standard arrangement of a pressure 



generator can be used with a flow diagnostic system 
that is adapted to provide the diagnostics outputs. 
The flow diagnostic system adapts itself to the 
characteristics of the pressure generator during the 
training mode and has a standard of comparison stored 
during the training mode that is available for 
comparison during the monitoring or operational mode. 
The standard of comparison can be adjusted as needed 
by a technician via the network. In each arrangement, 
the fluid flow meter provides a calibrated flow rate 
output and the diagnostic report indicates if the 
pressure generator is out of calibration. 

FIG. 8 is a flow chart of a process diagnosing 
the condition of a primary element. 

The condition of the primary element can include 
erosion or fouling of the primary element. The method 
or algorithm starts at 802. Digital pressure data is 
gathered in a training mode or time interval as shown 
at 804. A power spectrum of the sensor data, minus the 
moving average, is calculated at 806. The power 
spectrum obtained is identified as the training power 
spectrum at 808 and stored in non-volatile memory 810. 
After completion of training, the process moves on to 
monitoring or normal use . A further power spectrum of 
current pressure data, minus the moving average, is 
evaluated at 812, and the power spectrum so obtained 
in stored in memory 814, that can be either RAM or 
nonvolatile memory. At 816, the power spectrum Fi 
obtained during training is compared to the power 
spectrum Fi obtained during monitoring. If there is a 
significant difference between Fi and Fi which is 
indicative of a problem with the primary element, a 
primary element warning (PE Warning) is generated as 
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shown at 818. If the power spectrums Fi and Fi are 
sufficiently similar, then no primary element warning 
is generated. After the comparison at 816 and 
generation of a PE Warning, as needed, program flow 
moves to obtain new real time pressure data at 820 and 
the monitoring process moves on to a new evaluation at 
812, or the flow meter can shut down when there is a 
PE warning. The process in FIG. 8 can loop 
continuously in the monitoring mode to provide real 
time information concerning the condition of the 
primary element . 

In FIG. 9, a flow chart illustrates a process 
that provides diagnosis of both primary element (PE) 
and impulse lines (IL) . Program flow starts at 902. 
During a training mode illustrated at 904, pressure 
data, minus a moving average, is obtained and training 
power spectrum and training statistics are stored in 
nonvolatile memory as explained above. Next, impulse 
line diagnostics (such as those explained in Fig. 7) 
are performed at step 906 in FIG. 9. In FIG. 9, after 
impulse line diagnostics are performed, current 
impulse line statistics are compared to historical 
(training) impulse line statistics (as detailed in 
FIG. 7) at 908. If the comparison indicates a problem 
with plugging of impulse lines, then an impulse line 
warning is generated as shown at 910. If no problem 
with the impulse lines is apparent, then program flow 
moves on to primary element (PE) diagnostics at 912. 

At process 912, power spectral density for the 
current real time data is calculated (as explained 
above in connection with FIG. 8) . The current power 
spectral density is compared to the historical power 
spectral density at 914, and if there is a difference 



large enough to indicate a problem with the primary 
element, then a PE Warning is generated as shown at 
916 . If the differences in the power spectral 
densities are small, then no PE warning is generated 
as shown at 918. Program flow continues on at 920 to 
repeat the IL and PE diagnostics, or the flow meter 
can be shut down if there is a PE or IL warning until 
maintenance is performed. 

Any of the methods can be stored on a computer- 
readable medium as a plurality of sequences of 
instructions, the plurality of sequences of 
instructions including sequences that, when executed 
by a diagnostic application in a flow diagnostic 
system, cause the flow diagnostic system to perform a 
diagnostic method relative to a primary element and 
impulse lines coupled to the transmitter. 

The flow diagnostics system can also be used with 
a transmitter (not illustrated) which includes remote 
seals connected by flexible capillary tubes that are 
filled with a controlled quantity of isolation fluid 
such as silicone oil. The isolator arrangement permits 
the transmitter electronics to be spaced away from 
extremely hot process fluids which contact the remote 
seals. The flow diagnostics system can also be used to 
detect leaking and pinching off of capillary tubes 
using the diagnostic techniques described above to 
provide diagnostic data. 

The flow diagnostics system can also be used with 
a transmitter (not illustrated) which connects to taps 
near the bottom and top of a tank. The transmitter 
provides an output that represents a time integral of 
flow in and out of the tank. The transmitter includes 
circuitry, or alternatively software, that measures 
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the differential pressure between the taps and 
computes the integrated flow as a function of the 
sensed differential pressure and a formula stored in 
the transmitter relating the sensed pressure to the 
5 quantity of fluid in the tank. This formula is 
typically called a strapping function and the quantity 
of fluid which has flowed into or out of the tank can 
be integrated as either volumetric or mass flow, 
depending on the strapping function stored in 
10 transmitter. The transmitter can be located either 
near the bottom or the top of tank, with a tube going 
to the other end of the tank, often called a "leg." 
This leg can be either a wet leg filled with the fluid 
K in the tank, or a dry leg filled with gas. Remote 

j?y 15 seals can also be used with such a transmitter. 

K In one embodiment, the flow diagnostic system 

In includes a signal preprocessor algorithm that isolates 

I y 

3 signal components in the pressure data such as 

^ frequencies, amplitudes or signal characteristics 

p 2 0 which are related to a plugged impulse line or 

degraded primary element. The signal preprocessor 
provides an output to a signal evaluator in the 
difference algorithm. The signal preprocessor isolates 
a portion of the signal by filtering, performing a 
25 wavelet transform, performing a Fourier transform, use 
of a neural network, statistical analysis, or other 
signal evaluation techniques. Such preprocessing is 
preferably implemented in the difference algorithm or 
in a specialized digital signal processor. 
30 The signal components are isolated through signal 

processing techniques in which only desired 
frequencies or other signal characteristics such as 
amplitude are identified and an indication of their 
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identif ication is provided. Depending upon the 
strength signals to be detected and their frequency, 
signal preprocessor can comprise a filter, for example 
a band pass filter, to generate the isolated signal 
output. For more sensitive isolation, advanced signal 
processing techniques are utilized such as a Fast 
Fourier transform (FFT) to obtain the spectrum of the 
sensor signal. In one preferred embodiment, the signal 
preprocessor comprises a wavelet processor which 
performs a wavelet analysis on the sensor signal as 
shown in FIGS. 10, 11 and 12 using a discrete wavelet 
transform. Wavelet analysis is well suited for 
analyzing signals which have transients or other non- 
stationary characteristics in the time domain. In 
contrast to Fourier transforms, wavelet analysis 
retains information in the time domain, i.e., when the 
event occurred. 

Wavelet analysis is a technique for transforming 
a time domain signal into the frequency domain which, 
like a Fourier transformation, allows the frequency 
components to be identified. However, unlike a Fourier 
transformation, in a wavelet transformation the output 
includes information related to time. This may be 
expressed in the form of a three dimensional graph 
with time shown on one axis, frequency on a second 
axis and signal amplitude on a third axis. A 
discussion of wavelet analysis is given in On-Line 
Tool Condition Monitoring System With Wavelet Fuzzy 
Neural Network , by L. Xiaoli et al . , 8 JOURNAL OF 
INTELLIGENT MANUFACTURING pgs . 271-276 (1997). In 
performing a continuous wavelet transformation, a 
portion of the sensor signal is windowed and convolved 
with a wavelet function. This convolution is performed 
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by superimposing the wavelet function at the beginning 
of a sample, multiplying the wavelet function with the 
signal and then integrating the result over the sample 
period. The result of the integration is scaled and 
provides the first value for continuous wavelet 
transform at time equals zero. This point may be then 
mapped onto a three dimensional plane. The wavelet 
function is then shifted right (forward in time) and 
the multiplication and integration steps are repeated 
to obtain another set of data points which are mapped 
onto the 3-D space. This process is repeated and the 
wavelet is moved (convolved) through the entire 
signal. The wavelet function is then scaled, which 
changes the frequency resolution of the 
transformation, and the above steps are repeated. 

FIG. 10 is a graph of amplitude versus frequency 
versus time of a process variable signal such as the 
pressure data reading from a pressure transmitter. 
Data from a wavelet transformation of digital pressure 
data is shown in FIG. 10. The data is graphed in three 
dimensions and forms a surface 170. As shown in the 
graph of FIG. 10, the digital pressure data includes a 
small signal peak 172 at about 1 kHz at time ti and 
another peak 174 at about 100 Hz at time t 2 . Through 
subsequent processing by the signal evaluator, surface 
170 or portions of surface 170 are evaluated to 
determine impulse piping or primary element 
degradation. 

The continuous wavelet transformation described 
above requires extensive computations. Therefore, in 
one embodiment, a discrete wavelet transform (DWT) 
which is well suited for real time implementation in a 
control system or personal computer is performed. One 
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efficient discrete wavelet transform uses the Mallat 
algorithm which is a two channel sub-band coder. The 
Mallet algorithm provides a series of separated or 
decomposed signals which are representative of 
5 individual frequency components of the original 
signal . 

FIG. 11 is a block diagram of a discrete wavelet 
transformation. FIG. 11 illustrates^ an example in 
which an original set of digital/ pressure data or 
10 signal S is decomposed using a ^sub-band coder of a 
Mallet algorithm. The signal S yhas a frequency range 
from 0 to a maximum of fmx.f The signal is passed 
simultaneously through a first high pass filter 250 

i 

having a frequency range from 1/2 f max to f max / and a 

/ 

15 low pass filter 252 having /a frequency range from 0 to 
1/2 f max • This process is called decomposition. The 
output from the high pass filter provides "level 1" 
discrete wavelet transform coefficients 254. The level 
1 coefficients 254 , represent the amplitude as a 

20 function of time of that portion of the input signal 
which is between 1/2 f max and f max • The output from the 
0 - 1/2 fmax low pass filter 252 is passed through 
subsequent high pass (1/4 f max - 1/2 fmax) filter 256 
and low pass (0 - 1/4 fmax) filter 258, as desired, to 

25 provide additional levels (beyond "level 1") of 
discrete wavelet transform coefficients. The outputs 
from each low pass filter can be subjected to further 
decompositions offering additional levels of discrete 
wavelet transformation coefficients as desired. This 

30 process continues until the desired resolution is 
achieved; or the number of remaining data samples after 
a decomposition yields no additional information. The 
resolut/ion of the wavelet transform is chosen to be 
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approximately the same as the sensor or the same as 

the minimum signal resolution required to monitor the 

signal. Each level of / DWT coefficients is 

representative of signal ^amplitude as a function of 

time for a given frequency range. Coefficients for 

/ 

each frequency range are concatenated to form a graph 
such as that shown in/FIG. 10. 



In some embodiments, padding is added to the 
signal by adding data to the sensor signal near the 
borders of windows used in the wavelet analysis. This 
padding reduces distortions in the frequency domain 
output. This technique can be used with a continuous 
wavelet transform or a discrete wavelet transform. 
"Padding" is defined as appending extra data on either 
side of the current active data window, for example, 
extra data points are added which extend 25% of the 
current window beyond either window edge. In one 
embodiment, the padding is generated by repeating a 
portion of the data in the current window so that the 
added data "pads" the existing signal on either side. 
The entire data set is then fit to a quadratic 
equation which is used to extrapolate the signal 25% 
beyond the active data window. 

FIG. 12 is a graph showing signals output from a 
discrete wavelet transformation (DWT) such as that 
shown in FIG. 11. FIG. 12 illustrates an example 
showing a signal S generated by a pressure transmitter 
and the resultant approximation signals yielded in 
seven decomposition levels labeled level 1 through 
level 7. In this example, signal level 7 is 
representative of the lowest frequency DWT coefficient 
which can be generated. Any further decomposition 
yields noise. All levels, or only those levels which 




relate impulse piping or primary element degradation 
are provided. 

The algorithm evaluates the isolated signal 
received from the signal preprocessing and in one 
embodiment, monitors an amplitude of a certain 
frequency or range of frequencies identified and 
provides a diagnostic output if a threshold is 
exceeded. Signal evaluator can also comprise more 
advanced decision making algorithms such as fuzzy 
logic, neural networks, expert systems, rule based 
systems, etc . Commonly assigned U.S. Patent 
Application No. 08/623,569 describes various decision 
making systems which can be implemented in a signal 
evaluator and is incorporated herein by reference. 

The diagnostic application performs diagnostics 
related to the impulse piping or primary element using 
information derived from a differential pressure 
sensor. The following describes a number of 
embodiments for realizing a diagnostic application. 
The diagnostic application can provide a residual 
lifetime estimate, an indication of a failure, an 
indication of an impending failure or a calibration 
output which is used to correct for errors in the 
sensed process variable. 

A. POLYNOMIAL CURVEFIT 

In one embodiment of the present invention 
empirical models or polynomial curve-fitting are used 
to detect line plugging or primary element 
degradation. A polynomial -like equation which has a 
combination of input signals such as various 
statistical parameters can be used to detect primary 
element degradation or impulse line plugging. 
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Constants for the equations can be stored in a memory 
or received over the network. 
B. NEURAL NETWORKS 

The signal can be analyzed using a neural 
network. One such neural network is a multi-layer 
neural network. Although a number of training 
algorithms can be used to develop a neural network 
model for different goals. One embodiment includes the 
known Backpropagation Network (BPN) to develop neural 
network modules which will capture the nonlinear 
relationship among a set of input and output (s) . 

FIG. 13 is a diagram showing a simplified neural 
network. FIG. 13 illustrates a typical topology of a 
three- layer neural network architecture that can be 
implemented in the diagnostic application. The first 
layer 270, usually referred to as the . input buffer, 
receives the information, and feeds them into the 
inner layers. The second layer, in a three-layer 
network, commonly known as a hidden layer 272, 
receives the information from the input layer, 
modified by the weights on the connections 274 and 
propagates this information forward. This is 
illustrated in the hidden layer 272 which is used to 
characterize the nonlinear properties of the system 
analyzed. The last layer is the output layer 276 where 
the calculated outputs (estimations) are presented to 
the environment . 

FIG. 14 is a diagram showing a neural network 280 
used to provide a residual lifetime estimate. FIG. 14 
illustrates a schematic of a neural network which 
provides a residual life estimate 282 for a primary 
element or impulse pipe based upon a sensor 
(transmitter) signal 284. The sensor signal can be 



either a raw sensor signal or a sensor signal which 
has been processed through signal processing 
techniques . 

FIG. 15 is a graph 290 of residual life 292 
versus time 294. FIG. 15 illustrates a graph of 
residual life versus time and shows that an alarm 
level 296 can be set prior to an estimated failure 
time 298. This allows the system to provide an alarm 
output prior to actual failure of the device. 

C. THRESHOLD CIRCUITRY 

This embodiment uses a set of if-then rules to 
reach a conclusion on the status of the impulse piping 
or primary element. This embodiment may be implemented 
easily in analog or digital circuitry. For example, 
with a simple rule, if the signal drops a certain 
amount below a historical mean, an output can be 
provided which indicates that an impulse line is 
plugged or is in the process of becoming plugged. Of 
course, more complex rules can be used which use 
multiple statistical parameters or signal components 
of the sensor signal to provide more accurate or 
different information. 

D . WAVELETS 

With this embodiment , one or more of the 
decomposition signal (s) in a wavelet analysis directly 
relate to line plugging and are used to diagnose the 
transmitter . 

FIG. 16 illustrates a diagnostic fluid flow meter 
that has a pitot tube for a primary element . The 
primary element can comprise a simple pitot tube or an 
averaging pitot tube as illustrated. The pitot tube 
190 can be inserted through a tap 192 on a pipe 194. 
An instrument manifold 196 can couple between the 
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pressure generator (including pitot inlets and impulse 
pipes internal to the pitot tube 190) and a pressure 
transmitter 198. 

FIG. 17 illustrates a computer platform 1 that 
connects via an interface 2 to one of several Hart 
interchange Devices 4. Interface 2 can be an RS232- 
RS485 converter, an ethernet connection, an intranet 
or internet connection, or other suitable interface 
that communicates information to the computer platform 
1. The computer platform 1 is typically a personal 
computer located in a maintenance shop area that 
includes application software such as an Asset 
Management Solutions (AMS) software application from 
Rosemount Inc. Each Hart interchange devices 4 
couples to one or more pressure transmitters 6 via a 
termination panel 8 . The Hart interchange devices 4 
are coupled via a DIN rail or bus 10 to a control 
system 12 . A diagnostic application 14 as described 
above in connection with FIGS. 1-16 also resides on 
computer platform 1. Computer platform 1 provides a 
diagnostic report as explained above. The arrangement 
illustrated in FIG. 17 allows substantially all of the 
diagnostic software to run on computer platform 1 
(which is a small control system) rather than place 
additional overhead on control system 12 . The term 
"control system" as used in this application includes 
control systems such as control system 112 in FIG. 1 
which provide electrical feedback to a fluid 
processing plant as well as computers that perform a 
monitoring function such as computer platform 1, where 
the feedback to the fluid processing plant comprises 
human intervention based on a diagnostic report 
generated by the computer platform 1. 



Although the present invention has been 
described with reference to preferred embodiments, 
workers skilled in the art will recognize that changes 
can be made in form and detail without departing from 
the spirit and scope of the invention. For example, 
various function blocks of the invention have been 
described in terms of algorithms, however, many 
function blocks may be implemented in other forms such 
as digital and analog circuits, software and their 
hybrids. When implemented in software, a processor 
performs the functions and the signals comprise 
digital values on which the software operates. A 
general purpose processor programmed with instructions 
that cause the processor to perform the desired 
process elements, application specific hardware 
components that contain circuit wired to perform the 
desired elements and any combination of programming a 
general purpose processor and hardware components can 
be used. Deterministic or fuzzy logic techniques can 
be used as needed to make decisions in the circuitry 
or software. Because of the nature of complex digital 
circuitry, circuit elements may not be partitioned 
into separate blocks as shown, but components used for 
various functional blocks can be intermingled and 
shared. Likewise with software, some instructions can 
be shared as part of several functions and be 
intermingled with unrelated instructions within the 
scope of the invention. 



